Universality Class of Fiber Bundle Model on Complex Networks 
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We investigate the failure characteristics of complex networks within the framework of the fiber 
bundle model subject to the local load sharing rule in which the load of the broken fiber is transferred 
only to its neighbor fibers. Although the load sharing is strictly local, it is found that the critical 
behavior belongs to the universality class of the global load sharing where the load is transferred 
equally to all fibers in the system. From the numerical simulations and the analytical approach 
applied to the microscopic behavior, it is revealed that the emergence of a single dominant hub 
cluster of broken fibers causes the global load sharing effect in the failure process. 
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The fiber bundle model (FBM) has been studied for 
many years in order to explain a variety of failure phe- 
nomena caused by cascades pj. In the FBM, composed of 
N heterogeneous fibers put on a lattice, a fiber at the vth 
site is broken if the load a v is larger than the threshold 
value cr' ?l assigned to the fiber following a given proba- 
bility distribution function. When the fiber is broken the 
load which was supported by the broken fiber is shared 
among intact fibers following a load sharing rule. The 
two most frequently studied rules are global load sharing 
(GLS) 2], in which the load of a broken fiber is equally 
shared with all intact fibers in the whole system, and lo- 
cal load sharing (LLS) 0], which allows only the nearest 
intact fibers to carry the load of a broken fiber. 

Depending on which load sharing rule is used, the FBM 
has been shown to exhibit totally different behaviors: For 
GLS, the FBM has been found to have a phase transition 
toward the global failure as the external load per fiber (a) 
is increased beyond the nonzero critical point a c Q. It 
has been known that the avalanche size distribution for 
the GLS follows the power-law form with the universal 
exponent —5/2 0,0, and the universality class has been 
identified from the measurement of critical exponents • 
In contrast, the FBM under the LLS rule has been shown 
to belong to a completely different universality class; i.e., 
the critical value of the load approaches zero as N is 
increased following the form tr c ~ 1/ ln( N) || |j| . Re- 
cently, the transition between the GLS and LLS regime 
has been studied using modified load sharing rules [l(|. 

Until very recently, most studies of the FBM have been 
performed on regular lattice structures. In a general per- 
spective beyond the fracture of material, however, cas- 
cading failure triggered by overloading happens also in 
real- world network systems For instance, the 

recent blackout in the United States and Canada was 
caused by cascading breakdown of elements through a 
power grid. The scenario of major blackout is very simi- 
lar to the idea of the FBM; once one element fails, then 
the neighbor of the elements fail by the increased load 
from the failure. 

In this Letter, we numerically and analytically study 



the FBM subject to the LLS rule on various network 
structures: i.e., the Erdos-Renyi (ER) model of a ran- 
dom network the Watts-Strogatz (WS) model of a 
small- world network 01 > an d the static model of a scale- 
free network ^5| . Even though the model studied in this 
work obeys strictly the LLS rule, it is found that the 
FBM on complex networks exhibits completely different 
universality: the critical behavior, the avalanche size dis- 
tribution, and the form of failure probability function 
coincide with those of the FBM under the GLS rule. 

First, we briefly describe the FBM under the LLS rule 
on complex networks. A version of blackout scenarios for 
cascading breakdown of power plants from overloading is 
very intuitive to understand how the FBM on complex 
networks works. Fibers attached to the vertices of un- 
derlying networks act as power plants, and the external 
load on fibers can be regarded as the demand for elec- 
tric power. If the demand for electric power exceeds the 
capacity of a power plant, the power plant gets discon- 
nected from the network and the demand is transferred to 
neighboring power plants through the transmission lines, 
the edges of the network. In this model, the underlying 
network is rigid while the fibers attached to vertices are 
damaged. 

The local load transfer of broken fibers through the 
edges of the underlying network is governed by the LLS 
rule. Under a non-zero external load Na, the actual load 
a v of the intact fiber v is given by the sum of a and 
the transferred load from neighboring broken fibers. To 
systematically handle the local load transfer from broken 
fibers to intact fibers, we define the load concentration 
factor K v = <r v /a as K v = 1 + m j/kj, where the 
primed summation is over the cluster of broken fibers 
directly connected to v, rrij is the number of broken fibers 
in the cluster j, and kj is the number of intact fibers 
directly connected to j. A simple example is shown in 
Fig.^ which contains two clusters of broken fibers, mi — 
3, fci = 4, and m 2 = 2, fc 2 = 3. If < (1 + 3/4 + 
2/3)er, the fiber at v will be broken and join the clusters 
of broken fibers. We note that this is the generalization 
of the ID model (kj = 2) HI [HI 
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FIG. 1: The LLS rule applied for the FBM on a complex 
network. The vertices bound to broken fibers and intact fibers 
are denoted by the broken symbols and the empty circles, 
respectively. 



To explore all values of the external load Na, we in- 
crease a quasistatically starting from zero. In a finite- 
sized system, the infinitesimal increment S of a is ex- 



pressed as S — min„ ^ a with the minimization 

(min^) for all intact fibers. This condition is equivalent 
to the increase of a just enough to break only the weakest 
intact fiber, which is the minimal condition to trigger an 
avalanche. 

Our numerical scheme is as follows: The threshold 
value of the load ct*' 1 € [0,1] is assigned to each fiber fol- 
lowing the uniform distribution function [lij ] . Start from 
a = and repeat the following two steps: (i) Increase a 
by the infinitesimal increment 5. (ii) Following the LLS 
rule, break the fibers with er*' 1 < K v a iteratively until no 
more fibers break. For each increment of a, the size s(a) 
of the avalanche is defined as the number of broken fibers 
triggered by the increment. The surviving fraction x(a) 
of fibers is the ratio of the number of remaining intact 
fibers to N when the external load reaches Na, and thus 
is written as 



(i) 



We also measure the response function x, or the gener- 
alized susceptibility, by using 



dx 
da 



1 

7VA 



a<a<a+A 



(2) 



where we choose the small enough value A = 0.0005 for 
the numerical differentiation. The critical value a c of the 
external load, which is one of the key quantities of inter- 
est, is defined from the condition of the global breakdown 
x(a c ) = . 

Figure |21 displays the critical value a c and the suscep- 
tibility x f° r the FBM under the LLS rule on various 
network structures, such as the local regular network, 
the WS network, the ER network, and the scale-free net- 
works. Strikingly, we find that the critical behavior of the 
FBM on complex networks is completely different from 
that on a regular lattice. More specifically, while a c for 
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FIG. 2: (a) The system size (JV) dependence of critical points 
(<7 C ) for various networks with N = 2 s , 2 9 , . . . , 2 15 vertices, 
(b) The susceptibility for the networks with N = 2 14 . p 
and 7 are the rewiring probability in the WS networks and 
the exponent of degree distribution P(k) ~ fc -7 in the static 
model [15|. respectively. The data points are obtained from 
the averages over 10 4 (10 3 for N = 2 15 ) ensembles. 



the FBM on a local regular network vanishes in the ther- 
modynamic limit and is described by <r c ~ l/ln(./V) for 
finite-sized systems (see the curve for p = in Fig. |2 cor- 
responding to the WS network with the rewiring prob- 
ability p — 0), <t c for all networks except for the local 
regular one does not diminish but converges to a nonzero 
value as N is increased. Moreover, the susceptibility di- 
verges at the critical point following \ ~ (<?c — &)~ a '°, 
regardless of the networks, which is again in a sharp con- 
trast to the local regular network [see Fig. Hfb)]. The 
critical exponent 0.5 clearly indicates that the FBM un- 
der the LLS rule on complex networks belongs to the 
same universality class as that of the GLS regime al- 
though the load-sharing rule is strictly local. 

The evidences that the LLS model on complex net- 
works belongs to the universality class of the GLS model 
is also found in the avalanche size distribution P(s): 
Unanimously observed power-law behavior P(s) ~ s~ 5 / 2 
in Fig. |3 (a) for all networks except for the local regular 
one (the WS network with p = 0) is in perfect agree- 
ment with the same behavior for the GLS case 0. On 
the other hand, the LLS model for a regular lattice has 
been shown to exhibit completely different avalanche size 
distribution [1,0. 

The failure probability F(a) is defined as the proba- 
bility that the whole system is broken when an external 
load a, or less, is applied. In the LLS regime, the failure 



probability was studied to test the weak-link hypothe- 
sis [ijj Uji |2(j . In the test for the Weibull form, one can 
see in Fig. 13(b) that F(a)'s for complex networks fall on 
a common line which coincides with F(a) for the well- 
known GLS case, which is very much different for LLS 
on regular lattices 0,0, H3 [ see the inset in Fig.|3{b)]. 
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FIG. 3: (a) The avalanche size distribution P(s) and (b) 
the test of failure probability F for the Weibull form. The 
solid line in (a) indicate the result for the GLS regime having 
exponent —5/2. The inset in (b) shows the result for the 
LLS regime on the regular network (p = 0), which has clearly 
different form from others. 



Consequently, we again confirm that the LLS model on 
complex networks belongs to the same universality class 
as that of GLS. 

What makes the LLS model on a complex network have 
identical critical behaviors of the GLS model? In order 
to answer this question, it is helpful to investigate the 
microscopic details of the failure process. As the exter- 
nal load increases, the clusters of broken vertices form, 
grow, and merge into larger clusters, ultimately resulting 
in the emergence of a dominantly large cluster (we call 
it DLC henceforth). In case of a regular lattice, all small 
clusters of broken fibers are roughly in an equal condition 
because of the underlying regular (and thus spatially uni- 
form) topology. As the external load increases, all small 
clusters grow at roughly the same rate, and thus the DLC 
emerges abruptly. In contrast, the emergence of a DLC 
on complex networks is a gradual process because the 
DLC is formed at an early stage of loading and keeps 
growing continuously, as shown clearly in Fig.QJa). The 
above observation indicates that the growth of the DLC 
plays a dominant role in the failure of fibers on complex 
networks. 

Together with the early emergence and the gradual 
growth of the DLC, the small- world behavior [l^ in com- 
plex networks provides a reasonable qualitative expla- 
nation of the GLS-like behavior of the LLS model on 
complex networks. More precisely, the small- world ef- 
fect causes the average distance from the DLC to intact 
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FIG. 4: (a) The largest cluster size tul in the failure process 
versus the external load a. (b) The average distance d from 
the DLC versus the surviving fraction x. (c) The fraction 
n(x) of nearest intact fibers of the DLC as a function of the 
surviving fraction x. The symbols and the lines represent the 
numerical data and their fits to the curve of Eq. @, respec- 
tively. The data points are obtained for the networks with 
N = 2 14 and averaged over 10 4 ensembles. 



fibers to decay fast towards the value below two as the 
size of the DLC increases (or the surviving fraction x de- 
creases) [see Fig. 0Jb)], which implies that most intact 
fibers are very closely located to the DLC and thus the 
system behaves similarly to the GLS model with all fibers 
are separated by the unit distance to each other. 

Finally, we apply the mean-field theory to the growth 
of the DLC. Let us assume the situation that the system 
is at the ith step of the load transfer from the DLC. 
The surviving fraction is written as and thus the 
cluster is composed of N(l — x^) broken vertices (we 
assume that there exists only one cluster) and has k^ 1 
nearest neighbor intact vertices. Following the LLS rule, 
the load crW of the nearest intact vertices of the cluster 
satisfies the load conservation condition, which yields 



1 + - ajW 



-a, 



(3) 



where = few/JV. Assuming that the load threshold 
a th of a fiber is randomly distributed to the whole sys- 
tem following the cumulative distribution P cum ((T ), we 
establish the recursion equation for x, 



x (t+i) = x {t) _ n (t) 



1 - 



1 _ P (VrWl 
± 1 cum ) 

xW 



(4) 
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In order to solve the equations, we have to know the func- 
tional form of rS l \ which is difficult to determine analyt- 
ically. Instead, we assume the simple rational functional 
form with two fitting parameters a,b € [0, 1]. 



bx {t) (l 



,(*) 



)/(l-ax®), 



(5) 



which is motivated as a generalization of the GLS form 
(corresponding to a = b = 1). The numerical data of 
fit very well to the form © as shown in Fig. life). 

When the failure stops propagating, the fixed point 
of the dynamics described by Eq. (@J satisfies x^ t+1 ^ = 
x^ = x* . The uniform threshold distribution P cum (a) = 
a, together with Eqs. results in 



a(x*) 



bx*(l - x*) 
1 + (6 - a)x* 



(6) 



The critical value a c of the external load and the sur- 
viving fraction at a c are then easily obtained from the 
maximum of a(x*), yielding 



= bx*. 



-1 + yT+T 

b — a 



(7) 



Near the GLS regime, 6=1 and a = 1 — e, one obtains 



1/4 -e/8, z*~l/2- 



(8) 



which indicates that the solution for the GLS regime is 
revisited at e — in a good agreement with the numerical 
simulation where a c has always been found to be smaller 
than the GLS value a c = 1/4 (see Fig- EJ - In addition, 
combining Eqs. J2J and ©, we obtain the critical behav- 
ior of x* near a c , 



X 



dx* 
da 



-1/2 



(9) 



confirming the universal exponent —0.5 for the GLS 
regime. 

In conclusion, we have investigated the FBM on vari- 
ous complex networks including the WS network, the ER 
random network, and the scale-free network. From nu- 
merical simulations and analytical approach, it has been 
found that although the LLS rule is strictly applied, the 
critical behavior of failure exhibits the characteristics of 
GLS: The critical value of external load is finite in the 
thermodynamic limit, the divergence of the susceptibility 
is described by the same critical exponent as in GLS; the 
avalanche size distribution and the statistics of failure 
display the unique behavior of the GLS regime. 
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